A robust constraint on cosmic textures from the cosmic microwave background 
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Fluctuations in the cosmic microwave background (CMB) contain information which has been 
pivotal in establishing the current cosmological model. These data can also be used to test well- 
motivated additions to this model, such as cosmic textures. Textures are a type of topological defect 
that can be produced during a cosmological phase transition in the early universe, and which leave 
characteristic hot and cold spots in the CMB. We apply Bayesian methods to carry out a rigorous 
test of the texture hypothesis, using full-sky data from the Wilkinson Microwave Anisotropy Probe. 
We conclude that current data do not warrant augmenting the standard cosmological model with 
textures. We rule out at 95% confidence models that predict more than 6 detectable cosmic textures 
on the full sky. 
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Introduction. Precision measurements of 

anisotropies in the cosmic m.icrowave background 
(CMB) radiation have been instrumental in establishing 
the standard "ACDM" model of cosmology: that the 
universe is composed mostly of dark energy and dark 
matter, with structures seeded by nearly scale-invariant 
Gaussian density fluctuations. In addition to establish- 
ing ACDM, the CMB is also an ideal observable for 
determining if there are departures from this baseline 
model. 

In this paper, we present a novel algorithm for the 
analysis of CMB data from the Wilkinson Microwave 
Anisotropy Probe (WMAP) P] to search for the pres- 
ence of a class of topological defects known as cosmic 
textures [2]. Although textures (and other topological 
defects, such as cosmic strings) have been ruled out as the 
dominant source for the primordial perturbations [31 [3] , 
their production is inevitable in theories in which a non- 
Abelian global symmetry is broken [S] . Previous work [Hi- 
ll] presented evidence, based on the properties of a single 
feature in the CMB, that ACDM should be augmented by 
adding cosmic textures. Implementing Bayesian model 
selection using data on the full sky, we are able to put 
the texture hypothesis to a much more stringent test. In- 
corporating this extra information, we conclude that the 
WMAP 7-year data do not warrant augmenting ACDM 
with cosmic textures, and place constraints on theories 
giving rise to textures. Our algorithm is easily extendable 
to incorporate better data, multiple datasets, and a more 
complete theoretical understanding of the properties and 
evolution of cosmic textures. 
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Cosmic texture theory. The theory of cosmic tex- 
tures posits a phase transition in the early universe in 
which a non-Abelian global symmetry is broken. In 
an expanding universe, different regions of the universe 
can be out of causal contact, obstructing the symmetry- 
breaking phase transition from occurring in the same 
manner everywhere in space [5]. Therefore, a scale- 
invariant set of knots in the symmetry-breaking order pa- 
rameter inevitably form: these are cosmic textures. Sub- 
sequent to the phase transition, knots from the distribu- 
tion come into causal contact with their surroundings and 
undergo collapse [H [TDHT5] . Upon collapse, textures un- 
wind when the gradient energy of the field configuration 
exceeds the energy required to restore the global sym- 
metry. As the field re-orders, the energy of the texture 
configuration is released as an outgoing shell of scalar 
field radiation. 

The gravitational potential associated with a cosmic 
texture varies in time as it collapses and subsequently 
explodes. CMB photons passing through an evolving tex- 
ture will be redshifted if they pass through a collapsing 
texture, and blueshifted if they pass through an explod- 
ing texture 13J. Each texture unwinding event there- 
fore produces an additive hot or cold spot on the sky, 
which can be approximated as a disc whose angular size, 
Oc, depends on the distance to the texture unwinding 
event and whose amplitude, e = Stt^Gt]'^, depends on the 
scale of symmetry breaking rj. The temperature profile 
in the central region of an unwinding event situated at 
the Galactic North Pole can be approximated as [TS] 
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where 9 and (f) a-re Galactic co-latitude and longitude, 
respectively, and p = {0, 1}. Here p = corresponds to 



a hot spot and p = 1 to a cold spot. The form of the 
modulation for large is presently unknown; following 
Refs. O [7] , we match onto a Gaussian profile at the half- 
maximum radius, 0^, — yji/29c. 

The angular scale distribution is determined by the 
evolution of the cosmological horizon during the matter- 
dominated era |13| . which is fixed by the late-time cos- 
mological parameters of the ACDM model. In addition, 
each feature is equally likely to be: (i) hot or cold, and 
(ii) located at any point on the sky, allowing us to define 
the prior over the "local" template parameters as 
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where < ^o < "■, < '/'o < 27r, and we take 2° < 6'^ < 
50°. The lower limit on 6c results from the large power 
on degree scales in the CMB; the upper limit stems from 
the fact that templates with Oc > 50° are large enough 
to cover the whole sky and overlap themselves, rendering 
Eq.[T] invalid. 

Different theories giving rise to textures yield differ- 
ent predictions for the symmetry breaking scale and fre- 
quency of texture unwinding events; however, all mecha- 
nisms produce CMB modulations of the form described 
in Eq. [T] Observationally, theories giving rise to textures 
are therefore differentiated only by the expected number 
of detectable texture unwinding events on the CMB sky, 
Ns, and their amplitude, e. In our analysis, the back- 
ground CMB fluctuations dominate the definition of de- 
tectability. The prior probability Pr(iVs, e) is set by using 
simulations to determine the parameter space to which 
our algorithm is sensitive, as we will discuss shortly. 

The ACDM-|-texture model can therefore be fully de- 
scribed by: the standard ACDM parameters; a set 
of "global" texture parameters, mo = {Ns,e}, la- 
belling theories; a set of "local" parameters, nii — 
{p,0c,Oa,(l)o}i, describing each texture; and theoretical 
priors on these parameters, Pr(mo) and Pr(mi|mo). To 



test the ACDM-|-textures model against vanilla ACDM, 
we need only vary those parameters that are unique 
to the more complex model [M]. We therefore fix the 
ACDM parameters to their best-fit values from the anal- 
ysis of WMAP 7-year data (TS] (hereafter referred to 
as WMAP7). We will now describe the specifics of our 
search algorithm. 

Searching for textures. The fundamental question 
posed by this analysis is: are the WMAP7 data better 
described by the standard ACDM cosmological model or 
ACDM plus cosmic textures? The goal is to calculate the 
joint posterior distribution of Ng and e, given the avail- 
able data. Pure ACDM corresponds to Ng = 0. We avoid 
the a posteriori selection effects associated with postdict- 
ing an explanation for anomalous portions of the data 
(see Ref. [IS] for an in-depth discussion) by performing 
an analysis of the full dataset. This is important, given 
that previous evidence [7] for cosmic textures in the CMB 
was based on the analysis of a single anomalous feature, 
the so-called CMB Cold Spot [8]. 

Given an expected number of detectable textures over 
the whole sky, Ns, the actual number of detectable tex- 
tures, Ns , is drawn from a Poisson distribution with mean 
fskyNs, where /sky is the fraction of the sky covered by 
the observations. The full posterior probability distri- 
bution of the global parameters describing the texture 
model, e and Ns, is given by marginalizing the likelihood, 
Pr(d|mi, . . . mjv^ , e, Ns, /sky), weighted by the prior, over 
the (unknown) actual number of textures and their indi- 
vidual properties. This is an extremely challenging inte- 
gral to evaluate directly, but a good approximation to it 
can be found by identifying the regions of this parame- 
ter space in which the likelihood is appreciable and only 
including these contributions [T71[TH]. Extending this for- 
malism to also incorporate the global parameter e allows 
us to self-consistently combine the evidence that each 
candidate is a texture into a global constraint on the tex- 
ture theory. The resultant expression (cf. Refs. [T71[T5]) 
is 
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where Pr(e, TVs) is the prior (the properties of which are discussed below), Pr(d|iVs — 0, /sky) is the likelihood for 
ACDM (i.e. the likelihood assuming no textures), and Ni, denotes the number of regions on the sky, or "blobs," 
containing candidate signatures, each labeled by bi. The actual number of detectable textures Ns lies between 
and Ni). The quantity A''i^2-''"s is one when all indices take distinct values and zero otherwise: it generates all 
permutations of Ng textures located in Ni, blobs, assuming no more than one texture per blob. Finally, the quantity 
Pb- (e), defined as 
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is a patch-based evidence ratio evaluated in each blob: this is a measure of how much better ACDM plus a sin- 



gle texture fits the data than pure ACDM, considering 
only the data in blob bi. The factor of /sky appearing 
in the denominator accounts for the fact that we are re- 
stricted to detecting textures outside the sky cut. Unless 
the data provide strong support for the presence of a tex- 
ture, the evidence ratio penalizes this more complicated 
model through the larger volume of parameter space that 
must be considered in constructing the priors, thus self- 
consistently implementing Occam's razor. 
The likelihood for blob 6,- is 



Pr(d6jp, 0„ 00, 00, e, A^. = 1, /sky) = 
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where A'pix^b; is the total number of pixels in the blob, db. 
are the data points in the blob, and C;,. is the pixel-pixel 
covariance matrix using only pixels contained in the blob, 
which includes the fluctuations due to ACDM as well as 
instrumental noise and the effects of the beam. 

Locating texture candidates. To evaluate Eq. [3J 
we must first identify the most promising candidates in 
the map. We do so by employing the suite of spherical 
needlet transforms [TOtE] defined in Ref. [T7|. Filtering 
CMB temperature maps with spherical needlets yields 
information about both the position and angular size of 
interesting features. The statistics of the filtered field 
(established using 3000 simulated Gaussian CMB real- 
izations) can then be used to assess the significance of 
a candidate. Applying the needlet transform to texture 
templates (Eq. [I]) of various sizes yields a lookup table 
specifying the needlet whose response is maximal at each 
texture size. This table can then be used to identify 
peaks in a filtered input map with a texture candidate 
of a certain size. To minimize the number of false detec- 
tions, while not discarding potentially interesting signals, 
we determine a set of size-dependent thresholds (identical 
to those in Ref. [17]) using an end-to-end simulation of 
the WMAP experiment (see Refs. [531 [53]) containing a 
ACDM CMB as well as realistic foregrounds and system- 
atics that we cannot include in our likelihood function. 
The thresholds chosen restrict the number of candidate 
textures - by definition false detections - to be of order 
ten. All thresholds and parameters in the needlet trans- 
form are fixed at this point. 

Sensitivity testing. To determine our ability to de- 
tect textures given our thresholds, we generate a set of 
CMB maps from the WMAP7 best-fit power spectrum, 
and place textures of varying e and 9c in both a region 
with low and high instrumental noise (as the noise prop- 
erties of the WMAP experiment vary according to posi- 
tion on the sky). We find that, for 2° < 6*^ < 50°, the 
significance threshold is certainly exceeded (and there- 
fore a candidate identified) for e > 10"'*. For a favor- 
able realization of the background CMB and instrumen- 
tal noise, candidates are detected for e > 2.5 x 10~^ at 
scales 6c ^, 5° and e > 5 x 10~^ at somewhat smaller 



scales. We use e = 2.5 x 10~^ as a lower limit for de- 
tectable textures, and neglect the effect of 9c on our can- 
didate detection efficiency as it is far less important than 
the factor of 9'^ in Eq. ^ 

Calculating the texture posterior probability. 
Once the candidate textures have been identified, the 
posterior probability distribution Eq. [3]can be calculated 
by first evaluating the patch-based evidence ratio Eq. [4] 
for each blob using the MultiNest [Ml [25] nested sampling 
software. This requires calculating the inverse covari- 
ance matrix C^, , which is extremely memory-intensive 
at full WMAP resolution: the necessary storage capacity 
scales with size as 9c- We therefore employ an adaptive- 
resolution analysis, processing each blob at the highest 
resolution possible given its size and the available compu- 
tational resources. This removes the limitation on blob 
size of Refs. [T711IH]. 

The only remaining quantity to evaluate in Eq.[3]is the 
prior Pr(e, Ns). We choose a uniform prior for e between 
2.5 X lO-'^ < e < 1.0 x 10"". The lower bound is an esti- 
mate of what is detectable with our pipeline, determined 
by the simulations described above. The upper bound 
comes from requiring that the symmetry-breaking scale 
for textures, r/, is below the scale of cosmological infla- 
tion. To be consistent with the lack of observed B-mode 
polarization in the CMB [TS] US] , the scale of inflation 
must be less than approximately 10*^ GeV, constraining 
e to be less than roughly lO"** (this agrees with the prior 
of Ref. [6^). 

We adopt a uniform prior on Ng between < iVs < 10. 
The comoving density of textures produced in a phase 
transition depends on the particular texture model in 
question, and can be determined from simulations. The 
total number of unwinding events is obtained by inte- 
grating this density over the four-volume swept out by 
the CMB photons. For example, simulations [6] of SU(2) 
textures indicate that we can expect to have causal ac- 
cess to roughly 7 textures with 9c > 2° in the CMB. The 
number of these unwinding events which are then de- 
tectable is mainly a function of our particular realization 
of the background CMB. Our choice of a uniform prior 
accounts for our ignorance of both the precise theory giv- 
ing rise to textures and the precise number of detectable 
textures in the context of a specific theory. This allows 
us to compare ACDM and all models that give rise to 
textures. Under the assumption of a uniform prior the 
posterior is simply proportional to the likelihood; results 
for a different prior on Ng could be obtained easily by 
reweighting the current posterior. 

The significance required to favor the ACDM-|-textures 
model can be understood by evaluating Eq. [3] using a set 
of simulated evidence ratios Pbi(e), assuming that two 
candidate textures have been located in the data. The 
evidence ratios are chosen to be either low-amplitude 
and fiat in e (the case where each blob yields no sup- 
port for the texture model), or Gaussian with varying 
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FIG. 1. Regions containing 68% (dark blue) and 95% (light 
blue) of the posterior probability distribution, Eq. ^ for the 
hypothetical situations described in the text. Each case con- 
tains two hypothetical texture candidates. Top row: both 
candidates are ACDM only, and have low-amplitude, flat ev- 
idence ratios, P6j(e). Middle row: one candidate is "texture- 
like", and has a Gaussian evidence ratio whose amplitude 
increases from left to right. Bottom row: both candidates are 
texture-like; again, their evidence ratio amplitudes increase 
from left to right. 



amplitude (indicating varying degrees of support for the 
texture model). In all cases, the Gaussians are chosen to 
peak at the same value, e = 5xl0~^, and have the same 
standard deviation, cr = 5 x 10^^. The amplitudes of 
the Gaussian peaks are selected so that / pb- (e) Pr(e) de 
is 1/20, 1, or 20. These values are indicative of weak, 
intermediate, and strong texture signals, respectively. 

The posteriors for all combinations of the simulated 
evidence ratios are shown in Fig. [l] When none of 
the candidate features support the texture hypothesis 
(top row), the posterior is exponentially decreasing in 
Ns- In this case, we would correctly conclude that 
pure ACDM is strongly favored, and no constraints on 
e could be extracted. When one or two blobs produce 
a peaked evidence ratio (central and bottom rows), it 
becomes possible to make a detection. As the ampli- 
tudes of the evidence ratios are increased (left to right), 
the posterior begins to bulge, before ultimately becom- 
ing peaked. We would correctly conclude that the data 
favor ACDM-f textures over pure ACDM if a peak in the 
posterior at iV^ 7^ was sufficiently higher than the value 
of the posterior at Ng = 0. Comparing the central and 
bottom rows of Fig. [T] a detection can be made either 
in the case where there is a single strong candidate, or 
the case where there is a number of moderately strong 
candidates (provided each pbi{e) is peaked in the same 
range of e) . For the one- and two-strong-candidate cases 
(the centre- and bottom-right plots in Fig. fl]), the peaks 
of the posterior (after marginalizing over e) are 10 and 
1500 times that of the value at Ng = 0, respectively. 



Calculating the posterior for the end-to-end simulation 
of the WMAP experiment yields the constraints shown as 
solid and dashed lines in Fig. [2] This posterior resembles 
the top row of Fig. [I] is peaked at N^ = 0, and is not 
significantly different from the input priors on the global 
texture parameters. We therefore correctly conclude that 
the end-to-end simulation does not contain textures. 

Results and conclusions from WMAP. We per- 
form our analysis on the foreground-subtracted 94 GHz 
W-band temperature map from the 7-year release of 
the WMAP experiment [2S] (prepared by subtracting a 
model of known astrophysical foregrounds, as described 
in Ref. |27]). The W band has the highest resolution of 
the five measured by WMAP, with a full-width at half 
maximum of 0.22°. To minimize the effects of residual 
foregrounds, we apply the KQ75 mask, which yields a 
sky coverage of /sky = 0.706. The candidate textures are 
the same as those identified in Ref. [17] , minus one which 
lies outside our prior on 9c- The features range in size 
from 2° to 17.25°, and we are able to process seven at 
full WMAP resolution, two at half WMAP resolution and 
the largest at a quarter WMAP resolution. Although the 
lower-resolution computation of the likelihood for these 
three largest features does result in reduced accuracy, the 
impact on the overall posterior is minimal as the prior for 
textures of such large size is very low (cf. Eq. [2]) . 

Evaluating Eq. |3] yields the posterior for cosmic tex- 
tures in the WMAP7 data shown in Fig. [2] as dark- and 
light-blue regions. The posterior is clearly peaked at 
Ns = 0, and we find the marginalized constraint on the 
expected number of detectable textures to be Ns < 5.9 
(at 95% confidence). We therefore conclude that the 
WMAP7 data do not warrant augmenting ACDM with 
textures. The marginalized constraint on the scale of 
symmetry breaking is found to be 2.6 x 10~^ < e < 
1.0 X 10-^ (at 95% confidence). 

While the posterior is peaked at Ns = 0, there is also 
a clear difference between the WMAP7 posterior and 
that of the end-to-end simulation (over-plotted in Fig. [2| . 
Comparing the WMAP7 posterior to the example plots 
in Fig.jT] our result is also consistent with a signal that is 
present, but too weak to provide a detection. The differ- 
ent shape of the posterior is determined almost entirely 
by two features, located at {I = 185°, & = —79°) and 
{I = 209°, 6 = —57°) in Galactic coordinates, the sec- 
ond of which is the Cold Spot [3121]. As in Ref. [H], we 
use information from the multiple frequency bands of the 
WMAP instrument to confirm that there is no detectable 
residual foreground contamination in these features. This 
strongly motivates an analysis with better data, as will 
soon be provided by the Planck satellite [53] , or a better 
candidate-location technique, such as one utilizing opti- 
mal filters [30] . There is also the possibility of including 
CMB polarization data, as textures would not induce a 
polarization signal, unlike the primary CMB perturba- 
tions [9]). All of these efforts are currently in progress. 
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FIG. 2. Regions containing 68% (dark blue) and 95% (light 
blue) of the posterior probability distribution, Eq. pi for the 
WMAP7 data, along with the corresponding contours of the 
68% (solid line) and 95% (dashed line) of the posterior proba- 
bility for the end-to-end simulation of the WMAP experiment, 
based on a pure ACDM model. 



These and other tests will lead to better constraints on - 
or, if a signal is present, a confirmation of - the texture 
hypothesis. 
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